solve non-homogenous d.e's using method of variation of parameters